Hardcore dimer aspects of the SU(2) Singlet wavefunction 
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We demonstrate that any SU(2) singlet wavefunction can be characterized by a set of Valence Bond occupa- 
tion numbers, testing dimer presence/vacancy on pairs of sites. This genuine quantum property of singlet states 
(i) shows that SU(2) singlets share some of the intuitive features of hardcore quantum dimers, (ii) gives rigorous 
basis for interesting albeit apparently ill-defined quantities introduced recently in the context of Quantum Mag- 
netism or Quantum Information to measure respectively spin correlations and bipartite entanglement and, (iii) 
suggests a scheme to define consistently a wide family of quantities analogous to high order spin correlation. 
This result is demonstrated in the framework of a general functional mapping between the Hilbert space gener- 
ated by an arbitrary number of spins and a set of algebraic functions found to be an efficient analytical tool for 
the description of quantum spins or qubits systems. 
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I. INTRODUCTION 



Coupling a large number of spin s = 1/2 to a total spin 
S = state (singlet) leads to a very rich class of wave- 
functions which are the central pieces of many physical prob- 
lems ranging from Condensed Matter to Quantum Informa- 
tion. Indeed, the singlet state plays a major role in Quan- 
tum Magnetism especially in the low temperature properties 
of SU(2) Heisenberg frustrated antiferromagnets : below the 
finite spin gap energy scale, the physics is dominated by rather 
exotic singlet ground states (valence bond crystal, plaquettes 
state, spin-liquid, . . . ) and often not so well understood low 
lying singlet excitations'. It also attracted considerable in- 
terest since Anderson suggested that the Resonating Valence 
Bond (singlet) state (RVB) is the relevant insulating parent 
state of high-T c superconductivity 2 . Interestingly, these gen- 
uine quantum objects naturally emerge in the different context 
of Quantum Information because they are somehow paradig- 
matic of entanglement. Indeed, the simplest singlet (S = 0) 
wavefunction, namely the dimer state (1/V2)(|TJ.) — II t)), is 
a maximally entangled state made out two spin 1/2. For large 
systems, singlet wavefunction can naturally lead to massive 
entanglement^. 

Obtaining a simple and tractable description of singlet 
states when the number of spins becomes large, is thus a cru- 
cial issue. In that respect, it has been shown since the early 
days of quantum mechanics that Valence Bond (VB) states 4 
are a much more convenient and elegant approach to this 
question than expansions of singlets as linear combinations 
of Sf eigenstates. Indeed, Hulthen showed 5 that any singlet 
state can be written as a (non unique) linear combination of 
VB states, i.e. arbitrary range coverings of the N sites with 
dimers : 
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where = (1/V2)(| Ulj) - I UW) and v is a 

dimer covering of the system, defined as a partition of the 
{1,...,N} ensemble into N/2 oriented couples of sites 
. . . , (iN/2,3N/2)} (see figure E). 
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FIG. 1 : Two VB states corresponding to the coverings T> and T>' . The 
overlap diagram obtained by the superimposition T> U T>' involves a 
set a closed loops. 



However, while emphasizing the role of an apparently in- 
tuitive source of spin correlation and entanglement (namely 
the "bond dimer" (1/V2)(|U) - lit))), this approach leads 
to a non-orthogonal and massively overcomplete 6 description 
of singlet states that, in turn, destroys any interpretation of 
dimers as hardcore objects 7 . Paradoxically enough, this con- 
ceptually problematic "dimer fuzziness" (caused by the intrin- 
sic non-orthogonality of these objects) does not lead to major 
inconsistencies when neglected. Let us select three signifi- 
cant examples, from the contexts of Quantum Magnetism and 
Quantum Information, where this apparently ill-defined no- 
tion of bond occupation appears either implicitly or explicitly 
and is shown to give yet interesting and consistent results : (i) 
The Liang-Doucot-Anderson (LDA) wavefunction 8 convinc- 
ingly reinforces the idea that spin correlations can be driven 
by individual VB strength tuning, and in the same spirit Sand- 
vik 9 obtained very consistent results for the VB length distri- 
bution while mentioning the potential definition problem, (ii) 
The Rokhsar-Kivelson quantum dimer model 10 that considers 
dimers as hardcore and orthogonal objects is widely accepted 
as a good starting point to the study of SU(2) highly frus- 
trated (or more generally gapped) Heisenberg antiferromag- 
(1) net a 11 ' 12 , (iii) Alet et alJ^- recently introduced a VB Entangle- 
ment Entropy numerically shown to capture all the features 
of the von Neumann Entanglement Entropy but defined as if 
dimers where hardcore and orthogonal objects. 

In this article, we use a general functional framework to 
describe any singlet wavefunctions that allows to demonstrate 
that a VB occupation number testing the presence/absence of a 



2 



SU(2) dimer on a given pair of spins can be defined and com- 
puted unambiguously as an intrinsic property of any SU(2) 
singlet state. The implications are threefold : (i) We show that 
SU(2) singlets actually share some of the intuitive hardcore- 
like features of quantum dimers, (ii) We give rigorous basis 
for relevant albeit apparently ill-defined quantities introduced 
recently in the context of Quantum Magnetism such as VB 
Length Distribution 9 or Quantum Information such as VB En- 
tanglement Entropy 13 and give a systematic scheme to define 
consistently a wide family of quantities analogous to high or- 
der spin correlation, (iii) The functional mapping provides a 
convenient framework for further studies of correlated and en- 
tangled quantum spin or qubits systems and can be extended 
to arbitrary values of spin s and total spin S. 

The paper is organized as follow: in the first section we 
introduce and demonstrate the validity of the functional map- 
ping. In the second part the VB occupation number is ex- 
plicitly defined and shown to be the parent quantity for VB 
Length Distribution 9 and VB Entanglement Entropy 13 -. The 
final section is devoted to the generalization of the VB occu- 
pation number concept to arbitrary order VB correlation func- 
tions and the extension to higher spin s and total spin S values 
of the functional mapping. 

IL FUNCTIONAL MAPPING 

We introduce on each site i of the iV-spin system a real (or 
complex) variable n and associate to any dimer covering T> 
the N- variable function <px> defined by : 

ipT>(xi,...,x N ) = Yl d (xi,Xj), (2) 

where d(xi, Xj) is a function of xi and Xj that will be deter- 
mined later. Each function <£>z>({xj}) of the mapping corre- 
sponds to a pure VB state \ipj>) defined by (|T). The key prop- 
erties of the VB states are their non-orthogonality and over- 
completeness, both being closely related. Let us us recall the 
overlap rule that allows to compute the overlap 15 between two 
VB configurations T> and V, 

Ov,v> = (<pv\f>v>) = Ev,v>.2 ni ~ N/ \ 0) 

rii being the number of closed loops in the overlap diagram 
obtained by superimposing both configurations T> and V 
and ex> T>> the sign coming from the relative orientations of 
dimers in T> and V (see figure [TJ. The overcompleteness 
of the basis implies its non-orthogonality and on the other 
hand, all the linear relations between VB states are implic- 
itly encoded into the overlap matrix O : Because any sin- 
glet state can be expressed as a linear combination of ar- 
bitrary range VB states, the rank of O is the singlet sub- 
space size 7V = Nl/((1 + N/2)\(N/2)\). Its size being 
Af = N\/{2 N / 2 {N/2)\), all the overcompleteness of the basis 
is described by the Af—Afo independent singular eigenvectors 
of O, Hence, in order to obtain a faithful functional repre- 
sentation, it is necessary and sufficient to determine d(xi,Xj) 
and define a scalar product (( ( fiv({xi})\\ipv ({xi}))) acting 



on these functions that mimics the usual scalar product for 
VB states : ({<pi>{{xi})\\<pv'({xi}))) = 0x>,x>'. 
Let us show that d(xi,Xj) can be chosen as 

d\Xi)Xj) — \/~2 ^')' 

and 

{{<Pv(.{Xi})\\<PT»({Xi}))) = WdXi^XtpwiiXi}). 

x i = ±l/2 

(5) 

Given the form (|2j, <pv({xi}) X <pw({xi}) involves bond 
terms and each variable Xi appears exactly twice. Hence, this 
term can be represented graphically by the set of even size 
closed loops of the overlap diagram (see figure [TJ. Loops be- 
ing disconnected, ((<£>x>({£i})||¥>:D'({£i}))) is just a product 
of the individual loop (I) traces contributions, 

(<M{si})n<Mfo}))>= n *» <® 

igt>vd' 

with, for a loop of size p, 

t\ ^ d(xi 1 , Xi 2 ^d(xi 2 , ) ■ ■ ■ d(xi p , X{ 1 ) . (7) 

x i= ±l/2 

Choosing (0J implies that only two sequences, {x^ — 
{±1/2, =Fl/2, ■ ■ ■ , Tl/2}, give a non vanishing contribution 
to t{. As a consequence, t; = 2.2~ p< ^/ 2 .£; withp(Z) the size 
(number of sites) of the loop and e\ = ±1 the contribution of 
the loop I to the global sign Ed.v — Ylievuv £ i- Because the 
total size of the loops is just the number of sites N, the overlap 
((^({xiMcpv'ilxi}))) = e v ,v'2 n '- N ^ 2 as expected. 
Using (|5j, we showed that the set of functions 

ipv(xi,...,x N ) = \\ (xi-xj), (8) 

(t,3')ex> 

have the same linear properties as the VB states ([TJ and thus 
provides an algebraic mapping of the singlet subspace 16 . 

We emphasize that this mapping is general as it establishes 
a one to one correspondence between any VB state ([TJ and its 
functional counterpart (|8j : any subclass of VB states maps to 
its corresponding subclass both defining isomorphic Hilbert 
spaces. In the next section, we show (i) the particular role of 
the subset made of bipartite VB states, (ii) that valence bond 
occupation number is unambiguously defined using bipartite 
VB states despite non-orthogonality and overcompleteness. 

III. VALENCE BOND OCCUPATION NUMBERS 

Bipartition. The overcompleteness of the full VB basis pre- 
vents to define consistently a VB occupation number. We now 
show that such a definition is possible by restraining expan- 
sions of singlet states to a subset of V, yet overcomplete and 
spanning all the singlet sector. Let us partition the TV-spin 
system into two A r /2-spin subsets denoted A and B 11 -. We 
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FIG. 2: (a) Arrangement of N spins around a circle, alternating A 
sites and B sites, (b) In the definition d 14b of bond occupation 
operator, the limit for Xk is taken to a = 1/2 for A sites and b ^ a 
for ,6 sites except for x% and a!j for which the inverse convention is 
used : a if i (or j) is part of and b if i (or j) is part of .4. 



consider the subset V b of V containing the A/ b = (TV/2)! bi- 
partite VB states : 



(9) 



V also provides an overcomplete description of the singlet 
subspace : let us arrange the N spins on a circle in such a way 
that A sites alternate with B sites (see figure[2](a)) and choose 
an arbitrary singlet state \t/j). As shown by Rumeref al. 4 , there 
exist at least one decomposition of |^) as a linear combination 
of non-crossing*^ VB states. A simple inspection at figure [2] 
(a) shows that these non-crossing states are also bipartite. The 
occurence of a non-bipartite dimer would lead to the existence 
of an odd sized arc, hence causing at least one crossing. As a 
consequence any singlet state \ip) can be written as, 



Wi = E Xv\<pv)- 



(10) 



x>ev b 



This decomposition is not unique due to J\f h — M independent 
linear relations that reflect the overcompleteness of V b . 

Bond Occupation. We denote by n^j\ the occupation num- 
ber of a bond (bipartite or not) in a bipartite VB state 
\<Pv) ® and define it as, 



1 if (i, j) belongs to T>, 
otherwise. 



(ID 



For any singlet state \ip) decomposed on V b as (jTDJ, we show 
that the bond occupation 



is independent of the decomposition ( fTOb used, hence provid- 
ing a consistent intrinsic definition of a bond occupation num- 
ber for any singlet state. This statement is valid if any relation 
of overcompleteness implies its bond occupation counterpart : 



^ Ax> \<fv) = 



x>ev b 



]T Apn (lj) (|^>)) = 0. (13) 
vev b 



To demonstrate this point, let us use the functional mapping 
presented in the previous section and introduce the operator 



T>uj) acting on the functions of the mapping, 

%i)/({a*}) = 2 N/4 I™ {d Xi -d Xj )f({x k }), (14) 

with a or (3 = -1/2 (resp. +1/2) if i or j e A (resp. B) 
(see figure |2](b)). One can easily check that "Duj) have the 



following properties : (i) Linearity : T>uj\{f+g) = T^uj) /+ 
T>uj\g, (ii) For any V S V b , "Dn^ip-p — 1 if S V and 
T^(i,i)fv = if 4- Using (JHJ, l.h.s. of equation (fT3l 
becomes, 



x>ev b 



(15) 



(p,q)EV 



and applying T)uj\ to (fT~5T > immediately leads to the r.h.s. of 
(fT3l which concludes the demonstration. 

At this point, it is important to mention that nuj\ is rem- 
iniscent of a bond correlation Cy = (— 4/3)S.£.Sj but is es- 
sentially different. Indeed, ,-\(|V ! i>)) — ( { Pr>\C'ij\ i PT>) f° r 
any pure VB state. But, in the generic case of a superposi- 
tion ( [Tol l, both notions decouple due to the off-diagonal terms 
induced by the non-orthogonality of VB states : 



(^\C ij \i,)=nl ij) (m+ J2 ^vXv'MCijlw)- (16) 



E 

v,V 
v=tv' 



Illustration. For N = 6 (then 7V = 5, TV = 15 and 7V b = 
6) we define the bipartite VB states taking A = {1,3,5} 
and B = {2,4,6}. The overcomplete bipartite VB basis 
is: |a) = [1,2] [3, 4] [5, 6], \(3) = [1,6] [3, 2] [5, 4], | 7 ) = 
[1,4] [5, 2] [3, 6], \S) = [1,2] [3, 6] [5, 4], \e) = [1,4] [3, 2] [5, 6] 
and |C) = [1, 6] [5, 2] [3, 4]. The sole relation of overcom- 
pleteness writes |a)+|/3)+|7)-|<$)-|£:)-|£) = 0. We can check 
that e.g. the occupation of bond (5, 2) on the singlet state 
3 1 a) + 2 |C) that can be directly evaluated to 3x0 + 2x1 = 2 
is indeed independent of the decomposition : 3|a) + 2\Q = 
5|a) + 2\(3) + 2| 7 ) - 2\S) - 2\e) leads to5x0 + 2x0 + 2x 
1 — 2x0 — 2x0 = 2. On the other hand, the expectation 
value of (—4/3)82-85 is 5 and not 4 = 2 2 for this state. 

VB Length Distribution & VB Entanglement Entropy. Us- 
ing this result any quantity built out from bond occupa- 
tion/vacancy scheme on a singlet state \tp) can be defined 
unambiguously. The VB Length Distribution computed in 9 
writes, 



vev b 



E*A(|V>» 



with £«(M) = E ^WM). 



with 
length I 



(17) 

where £i(\ip)) is the total occupation of bonds with length I. 
A definition of the VB Entanglement Entropy introduced in 13 
as an equivalent of von Neumann Entanglement Entropy mea- 
surement of the bipartite entanglement between subsystems fl 
and Cl spanning the all system is, 



°n,n 



(|V»=ln2x 



2WIV>». 



(18) 



crossing 
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IV. EXTENSIONS AND GENERALIZATIONS 

VB correlation functions. Comparing (0]i and (fl4l exhibit 
a duality in the mapping between the dimer representation 
(Xi — xj) and a dimer susceptibility (d Xi — d Xj ). It is then 
straightforward to introduce p— order VB occupation num- 
bers (or VB correlation functions), similar to 2p— spin cor- 
relation functions, by defining linear operators in the spirit 
of (fl4l > containing more than one bond susceptibility (d Xii — 
d Xjt ) . . . {d Xip — d Xjp ). The analogous of equation ( TToT l for 
such an operator, shows that both quantities are equivalent for 
a pure VB state but becomes distinct yet consistently defined 
for a VB superposition. While much simpler to compute, we 
infer that VB correlation functions provide the same insights 
than usual correlations functions when computed on a given 
singlet state. 

Higher total spin sectors for s — 1/2. Following Hulthen's 
scheme 5 higher total spin S sectors with fixed total S z can 
be (overcompletely) represented by introducing in the cover- 
ings (Q3 an arbitrary range totally symmetrized state of 25* 
spins {ii, . . . ,«2s}s,- The functional mapping can be ex- 
tended to these case in the spirit of using the generating 
function G s {n, {x ip }) = fj^i (l + V x i P ) and representing 
{k, ■■ ■ ,«2s}s, with 

asA x *i>---> x *s) ~ y (2S)\ (S+S z )\ V " y ) M=0 ' 

(19) 

In particular X2) = {l/\2)[xi + Xj) represents the 

S z = bond triplet. Mixed total spin state such as the 
Neel state can be represented as well in this framework : 
If we denote A the 7V/2-site sublattice with up spins, it is 
easy to check that the corresponding Neel state is represented 

b y u a = Yl lt zA x i since ' usin s ©' ((^({xiDHru)) = 

(X/ 2 )JV/4 for any bi p artit e £>. 

Singlet sector for higher spins s. As noticed by Tasaki 19 , 
VB concept can be extended to higher values of s : [i,j) s — 
73=71 ££*=-. I - 8» +a z ). The overlap Ov,v © 

becomes {^vWv) = £x>,x>'-(2s + \) ni ~ N / 2 with e-p,v = 



1 for integer s. It can be readily checked that choosing 
d a (x i} Xj) = -j£-f(-i)"- x '6 Xii - Xj and enlarging the sup- 
port of the sum defining the scalar product © to Yl x r =-s 
leads to ({<pv({xi})\\<pw({xi}))) = (fv\fv) thus provid- 
ing the extension of the mapping. The proof that s-spins VB 
states provide an overcomplete description of the singlet sec- 
tor (if N is large enough to ensure that the number of cov- 
erings is larger than size of the singlet sector) is beyond the 
scope of this article^. 



V. CONCLUSIONS 

Any SU(2) singlet state can be characterized by a set of 
bond occupation numbers : despite the fuzziness induced by 
the overcompleteness of the VB basis, singlet wavefunctions 
conserve intuitive hardcore dimer features such as the pos- 
sibility of deciding whether a SU(2) dimer live on a bond 
or not. This genuine quantum property formalizes the sug- 
gestive idea behind the LDA wavefunction 8 that VB bond 
strength can drive the correlations properties of the singlet 
wavefunction and gives insights to the question of the inter- 
play between non-orthogonality, overcompleteness and hard- 
core properties. Moreover, this allows consistent definitions 
of quantities relevant for Quantum Magnetism (VB Length 
Distribution 9 and VB correlations 14 of any order) or Quantum 
Information (VB Entanglement Entropy 13 ) emphasizing the 
intuitive role of dimers as elementary bricks to build correla- 
tions or entanglement. We note that all these quantities are 
particularly suitable for numerical calculation in the frame- 
work introduced by Sandvik 9 . From a methodological point 
of view, the functional mapping used to derive this result ap- 
pears as a suggestive tool to describe quantum spin or qubits 
systems and can be extended naturally to higher spins s and 
higher total spin S sectors. 
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tionale de la Recherche (France). 
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